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■r^ ' Abstract 

o: 

^D ' We address the question of radiation emission from a collapsing star. We consider the simple model of a spherical 

vH I star consisting of pressure-free dust and we derive the emission spectrum via a systematic asymptotic expansion of 
f5 ■ the complete Bogolubov amplitude. Inconsistencies in derivations of the black body spectrum are pointed out. 

^-5 ; 

00 ■ I. INTRODUCTION 

^ . More than a quarter of a century has passed since Hawking's remarkable suggestion that a star collapsing to a 

■^ ■ black hole gives rise to radiation emission at a steady rate characterized by the black body spectrum (Hawking 1975). 

t^^ That particle creation takes place is not in itself surprising. Let us take the case of a spherical star with radius 

^^ initially larger than its Schwarzschild radius that eventually collapses contracting to a point (according to classical 

gravitation). If we consider a quantized scalar photon field the quantum spaces of the in states (before the initiation 

of collapse) and of the out states (after collapse is completed) are certainly different. Hence particle creation clearly 

^~j [ takes place and is determined by the Bogolubov a{uj, u') and (3{uj, u') amplitudes. The fact that the spectrum is that 

— «. of a black body is indeed noticeable, and ties up with the somewhat earlier results on black hole thermodynamics. 

^^ From the mathematical point of view the black body result is due to the special behaviour of the photon modes near 

I , the stellar surface just before the horizon is formed. Hawking in his derivation makes heavy use of the asymptotic 

rTp.' (i.e. near the horizon) form of the modes. The singular behaviour of the modes in this regime has in turn given rise 

. . I to various statements in the literature that are not strictly correct. For example there have been references to "parts 

^ • of a{u),uj') and /3(w, w') that relate to the steady-state regime at late times" ((DeWitt 1975), p. 327). However the 

k/rH " Bogolubov amplitudes are global constructs and the distinction between early and late times does not make strict 

''jj . sense. Hawking too talks about particle production that depends on the details of the collapse and that such particles 

d ' "will disperse", thus leaving only the thermal part at late times (Hawking (1975), p. 207). The point of view we are 

adopting in this paper is the one dictated by quantum-mechanical orthodoxy, namely that particle production is a 

global process ((Hawking 1975) p. 216) and that one should start with the full standard expression for P{uj,uj'). Of 

course at some point in the mathematical handling of the amplitude P{ll!, uj') the special role of the horizon will show 

up. The above remarks are certainly not meant to imply that one cannot calculate local quantities like (T^i/); such 

quantities certainly behave very differently during the various phases of the collapse. 

Shortly after Hawking's work on black hole evaporation two classic papers (Fulling and Davies 1976) and (Davies 
and Fulling 1977) were published (the latter shall be quoted as DF in what follows). The authors demonstrated 
an illuminating analogy, physical as well as mathematical, between gravitational collapse and the seemingly rather 
different problem of a perfect mirror starting from rest and accelerating for an infinite time. The renormalized matrix 
element of the r„„ component of the energy momentum tensor, defined as 

Tun = {du(t>f 

is calculated in DF. The result is that asymptotically for t —* oo 



IT^\ ^ (1) 



(see e.g. equation (4.16) of (Birrell and Davies 1982) where a comprehensive review of particle production from both 
mirrors and black holes is presented; k is a constant characterizing the trajectory.) Equation (1) shows that there is a 
constant energy flux at late times, analogous to the thermal energy flux found in (Hawking 1975). Davies and Fulling 
(1977) also calculated the Bogolubov amplitude following Hawking's lines and arrived at the black body spectrum 
relevant to an accelerating mirror 
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l^(-''^')l^w;:;w^ (2) 

POO 

n{uj)^ / dLo'\(3{uj,uj')f (3) 

Jo 

Certain technicalities in the latter calculation have been elucidated in a previous paper by the author (Calogeracos 
2002, hereafter referred to as I). In the present paper we intend to show that similar points may be raised in connection 
to the standard black hole literature (see (Birrell and Davies 1982), chapter 8 and references therein). There are steps 
in the derivation of the black body spectrum that both obscure the mathematics and create confusion regarding 
the physics of the problem. Our findings, summarized in the concluding section, in no way do they diminish the 
astounding intuition shown by the early workers on the subject. 

In section 2 we examine the simple model of a collapsing spherical star consisting of pressure-free dust and we 
review the fundamental features of the collapse. The problem is widely treated in the literature and some of the 
results are included so that the paper be self-contained. In section 3 we write down the photon modes and the full 
expression for the Bogolubov amplitude P{u},u}'). In section 4 we show that the asymptotic behaviour of the latter 
for large uj' is 

Piu;,uj') « {oj')-'' + O Uoj')-'') (N > 1) (4) 



The l/co' in (2) leads to a logarithmic divergence in (3) and this signals the production of particles at a finite rate 
for an infinite time. One often refers to the ultraviolet divergence by saying that large to' frequencies dominate. Our 
analysis emphasizes two points, stressed also in I: (a) the thermal result depends crucially on the behaviour of the 
photon modes near the horizon, (b) for a consistent derivation one must consider the whole collapse and not just the 
late phase. The truth of statement (a) is usually taken as common knowledge. However the importance of statement 
(b) is often not appreciated. Comparison of our approach with the somewhat more conventional one is presented in 
Appendix A. 

The ultraviolet divergence previously mentioned is widely taken to signify that one cannot fully resolve the problem 
within the context of classical gravitation. It is well known that the constants we have at our disposal are c — 3x 10^° 
cm-s~^, /h = 1.05 X 10"^"^ g-cm^-s""'^, G = 6.67 x 10~^ cm^-s~^-g~"'^ and that we can form three quantities with 

dimensions of mass, length, and time respectively: M = (/ic/G) = 2 x 10^"'' g, L ~ (^/iG/c^) = 1.6 x 10"^"^ 

cm, T = (/iG/c"'') = 5 X 10^'''* s. It is clear that when the radius of the contracting star becomes of order L then 
quantum gravitational effects have to come into play. Such matters are not discussed here. Nor do we touch upon 
the profound problem of loss of information during collapse (see e.g. Preskill 1992, Heifer 2003). 
Note: in what follows /i = G = c = 1. 

II. A COLLAPSING SPHERE OF DUST 

We consider a collapsing spherical star consisting of pressure- free dust and follow the approach of (Weinberg 1972), 
chapter 11, sections 8 and 9. Since each dust particle falls freely the spacetime geometry inside the star is most 
appropriately described in a comoving frame. Let M be the mass of the star. The metric reads 



ds^ ^dt^ -R^{t) 



dr 
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1 — kr^ 



+ r^d0^ + r^sm^0dip 



(5) 



In the comoving frame each dust particle is labelled by a unique set of (r, 6, ip) and the time t stands for the proper 
time registered by the particle in question. The radius of the star is specified in the comoving frame by r = a where a 
is by definition constant. The quantities k and R{t) refer to details of the model and their meaning shall be explained 
shortly. The collapse is initiated at time t = 0, and we normalize R(t) by requiring 



R{0) ^ 1 
Outside the star spacetime is described by the Schwarzschild metric 

ds^ =fl- ^\ dp -(l- ^) df^ - fHe^ - f2 sin2 9d^^ 

The match of the two metrics at the surface imphes (Weinberg 1972) 

f ^ rR{t),e ^e,ip^(p 



(6) 



(7) 



(8) 



Relation (6) together with the first of (8) imply that the quantity a stands for the radius of the star in Schwarzschild 
spacetime at the moment when collapse starts. Thus the radius of the star in Schwarzschild coordinates is given by 
R = aR. The constant k is proportional to the star's density and is related to the mass via 



2M = ka-^ 
At times we shall consider the case of an initially very large star a » M, which implies via (9) 

ka^ « 1 



(9) 



(10) 



The relation between t and i is a rather more intricate one and is treated in detail by Weinberg op. cit. (see the 
remark after (14) below). It has already been mentioned that since the pressure inside the star vanishes the motion 
of the surface corresponds to that of a free falling particle in the gravitational field of a body of mass M . The proper 
time t registered by an observer located on the surface of the collapsing star and the quantity R are best given 
parametrically in terms of the cycloidal variable rj (Weinberg 1972): 



R{r]) = aR{7]), R{t]) = - (1 + cos r;) 



(11) 



t{v) 



2Vk 



{rj + sin rj) 



(12) 



where < rj < n. Collapse starts at ry = and is completed at r] = tt. Thus according to (12) the total time of the 
collapse as measured by an observer comoving with the star's surface is given by tt/ ( 2\/fc ) • The black hole is formed 
when R = 2M. According to (11) this happens at a value rjo of the parameter given by 
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= cos 



aVk 



(13) 



The Schwarzschild time t for a particle on the stellar surface is given in terms of 77 by the expression (see (Chandrasekhar 
1983) chapter 3, section 19) : 



t = 2a\/\ - ko? 



1 



(j] + sin 77) + ka rj 



2M In 



tan^ 



tan ^ 



tan ^ - tan 2 



(14) 



Relations (12) and (14) parametrically connect the time t in the comoving frame to the Schwarzschild time t on the 
stellar surface. Note that as 77 ^ 770 from the left (i.e. before the black hole is formed) t diverges. This reflects the 
well-known fact that collapse appears to the Schwarzschild observer to take an infinitely long time. 
Radial null geodesies in Schwarzschild spacetime satisfy (see (7)) 



dp 



1 



(1 



2M^ 



rdf' 



(15) 



This leads to the definition of the Regge- Wheeler radial coordinate r* 

r*-f + 2Mln(^-l) (16) 

so that (15) may be written as dP — dr*^ ~ 0, and to the introduction of the incoming and outgoing Eddington- 
Finkelstein coordinates 



v^t + r*,u^t-r* (17) 

The quantities v and u label photon trajectories in Schwarzschild spacetime. 

In order not to confuse the reader with the proliferation of radial variables let us clarify: R stands for the stellar 
radius in Schwarzschild coordinates, R* is the corresponding Regge- Wheeler coordinate, R is the auxiliary function 
given by the second of (11) and appearing in (5), f is the radial Schwarzschild coordinate of an arbitrary point outside 
the star, r* the corresponding Regge- Wheeler coordinate, and r the radial coordinate inside the star in the comoving 
frame. 

A property that one must establish is that as the star approaches collapse 

J- - 1 ~ Ae-*/2Af (^8) 

2M ^ ' 

where A is a calculable constant. This expression is of course widely known. It agrees with the remark following (14) 
that as far as the Schwarzschild observer is concerned the collapse takes an infinitely long time. A similar exponential 
appears in the calculation of the red shift as observed by the Schwarzschild observer; see (Weinberg 1972), p. 348 
(the asymptotic form of the mirror trajectory given by DF, equation (2.3), has a form identical to that of (18)). 
To establish the validity of (18) and also calculate A one takes logarithms of both sides and uses (11), (14). The 
logarithmic divergences cancel and A is expressed in terms of r]o (see Appendix A). The quantity A depends on the 
parameters a and M of the collapse and does not turn up in the expression for the spectrum. Note that relations (17) 
define parametrically (via rj) a function 

u^fstiv) (19) 

and its inverse 

v=Pst{u) (20) 

which give the trajectory of the stellar surface in terms of u, v coordinates. 

In the limit rj ^ rjQ both R* and i are singular (as is obvious from (14) and (16)), however their combination in the 
Eddington-Finkelstein coordinate v is analytic. The cancellation of the logarithmic singularities is hardly surprising 
since vq must have a well-defined value labelling a null line; for details see Appendix A. For given values of M and a 
the value of 779 is readily calculable via (13) and then so is vq after some algebra (setting 77 = ?7o in (81)). Thus vq — v 
admits a Taylor expansion near ryo ^nd of course so does 2M — R which is analytic for all < ry < tt. Writing 

2M~R^Ci{rjo-r^),vo^v^C2{m-rj) (21) 

we can read the positive constants Ci, C2 off (83), (82). Then 

i?-2M~C(wo-w),C = Ci/C2 (22) 

Combining the above with (18) we get 

.c.z;o-^e-*/2M (23) 

For rj ^ rjo we can combine the two relations (17) and write 

t c -^ (24) 

Then (23) reads 

v,,vo-Be--/''\B^^ (25) 

This defines Pst{u) in (20). 

In what follows we shall need the metric (7) expressed in terms of Kruskal coordinates; see e.g. (Townsend 1997), 
(Misner, Thorne, Wheeler 1973). We introduce 

v/4M 

U = -£;e-«/4^, V = ——- (26) 

E 



and we transform the metric to 



ds^ = ^^e-^/^^dUdV - f^dn (27) 

r 

where r is given in terms oi U ,V implicitly via 

UV = ^_^e-72M (28) 

2M ^ ^ 

At early times spacetime is taken to be essentially flat (cf remark preceding (35)) and the metric is simply expressed 
in terms of advanced and retarded Eddington-Finkelstein coordinates 

ds^ = dudv - r'^dn (29) 

We now turn to ray-tracing and consider light rays obeying the condition that they be reflected at the centre of 
the star r = 0. The rationale for the boundary condition at r = is reviewed in the next section. An incident ray 
corresponding to a certain v at early times upon reflection becomes an outgoing ray corresponding to a certain u at 
late times, thus defining a function u = f(v). The inverse function is given by w = p(u). Note that the function 
u ~ f{v) is by construction quite distinct from the function u = fst{v) (19) (the former involves reflection of the 
ray at the centre of the star whereas the latter at the surface). The importance of ray tracing lies in the fact that 
once we determine the function f{v) (or its inverse) we can immediately construct the photon modes via (45), (46). 
Note also that these expressions involve the Eddington-Finkelstein coordinates rather than the Kruskal ones. It is 
convenient to exhibit collapse via a Kruskal diagram, where radial null lines are at 45° on the plane of the paper. The 
line 1(1) represents a ray that crosses the star's surface at point B, is reflected at the centre of the star at spacetime 
point C, and after reflection crosses the star's surface at A (reflected ray labelled by /'(!))■ By definition ray 1{1) 
is the last incoming ray that manages to escape before the black hole is formed, and corresponds to an incoming 
Eddington-Finkelstein coordinate vh- At A the star's surface crosses the Schwarzschild radius, and the coordinates 
of A correspond to the value rjo of the cycloidal parameter 77 in (11), (12). The extension of the line CA is the future 
horizon and corresponds to the value U — oi the Kruskal coordinate. Let l'{i) be an outgoing ray slightly preceding 
1(1), corresponding to an incoming ray l{i). We also trace the incoming light ray ?n(l) which becomes 1{1) upon 
refiection on the star's surface at A without entering the star (this is an auxiliary ray and the reflection mechanism is 
an entirely hypothetical one, not connected to partial reflection that may take place on the star's surface). Similarly 
let m{i) be an incoming ray slightly preceding m{l). The problem is to relate the u of the outgoing l'{i) to the v 
of the incoming l{i), thus determining the function u = f{v). We review the argument in (Hawking 1975); see also 
(Townsend 1997), p. 125. Recall ((Townscnd 1997), p. 29 for proof) that U is an affine parameter for the null geodesic 
U =constant. (This means that L • DL^ = rather than L • DL^ ocL^.) Instead oi U = — e-«/4M ^j^g j^g^y g^jgQ ^gg 
any parameter A — — _Ee~"'^*^, E being some constant; we exploit this freedom and will comment on the value of 
E later on. Let T be a null vector on U — parallel to the horizon and pointing along the radial spatial direction, 
and let N be a null vector pointing to the future and along the radial spatial direction and satisfying N • T = —1. 
The above requirements are consistent (see (Townsend 1997), p. 109) if N is parallely transported along the geodesic 
defined by T. Consider a neighbouring null geodesic labelled hy U ~ —e where is small. This corresponds to constant 
u according to (26): 

e = £;e-"/4A^ ^ In e = In £; - -^ (30) 

AM ^ ' 

The two geodesies are thus connected by the displacement vector — eN. We parallely transport the pair N, T to the 
point where the [/ = geodesic cuts the trajectory of centre of the star. We then reflect N and T as in figure 2 
and parallely transport the new pair to past infinity where the stellar radius satisfies (35). Then spacetime becomes 
essentially flat, the metric is given by (29) and the affine parameter is given by v for a general geodesic l{i) (figures 
1, 2) and by vh for the special geodesic /(I). We thus have 

£ — Vu ~ V (31) 

Comparing (30) and (31) we obtain 

vh-v = Ee-^l^^^ 

or 



v^VH- Se-"/4M (32) 

(the ~ symbol above refers to the fact that this approximation is vahd near the horizon). Note that the preceding 
argument involving N, T can be equally well applied to the incoming rays m{i), m{l) which upon reflection on the 
surface of the star coincide with l'{i), l'{l)- In particular the afhne distance between 1{1) and l{i) equals the distance 
between ?ti(1) and m(i). This implies that the quantity E appearing in (32) is equal to the B appearing in (25). (The 
identification of E with B does not affect the derivation of the black body spectrum and we shall keep E in notation.) 
Note also that the quantity 

d=VQ-VH (33) 

is calculable; see appendix A, the remark following (89). Inverting (32) we get 

We now turn to the early stage of the collapse (figure 3). We assume that the collapse is initiated at t ~ (and 
that for i < the star is held stable by some external means) . Let us assume large initial stellar radius 

a» M (35) 

Expanding (11), (12) for early times, i.e. small r], we obtain 

Vk 4 

Hence the radius as a function of time at the early stage is given by 

Recalling the definition (9) we observe that the star's surface initially contracts according to non-relativistic kine- 
matics with acceleration M/a^ (in accordance with Newton's law of gravity). Given the assumption (35) the initial 
gravitational field is sufficiently weak so the light rays in the Schwarzschild geometry are straight lines. Let us consider 
a.i — R diagram ( t, R being Schwarzschild coordinates; at early times when (10) is satisfied there is little difference 
between radial and Regge- Wheeler coordinates) and take an incoming light ray q which has u — 0. The reflected ray 
q' has u ~ Q and takes time a to travel from the centre of the star to the surface. During this time the gravitational 
field is still weak. This may be seen by setting t = a in (36); then (10) implies that by the time the ray emerges from 
the star's surface R still is almost equal to a. To summarize the function u — f{v) has the properties 

/(f) = f = 0, initially (37) 

f{v)^-4Mln[^^^^],v^VH (38) 



III. CONSTRUCTION OF THE IN AND OUT STATES 

We follow (Birrell and Davies 1982), (Brout et al. 1995). A scalar mode corresponding to angular momentum I 
and satisfying the Klein-Gordon equation is written in the form 

i^i = (^VI^ry\i{r)Yr{e,^) (39) 

where 



d^ d 



2 



dP dr*^ 



Vi{r)]<l)i^O (40) 



with 



ViM = |i-^V^ + ^^l («) 



r J \ r 

As far as conceptual purposes and technical details are concerned it suffices to restrict ourselves to s waves (following 
the references in the beginning of this section). Note that according to (41) there is a centrifugal barrier even for s 
waves which does not affect either the high frequency modes or the modes that are present at the early stages of the 
collapse (when the gravitational field is weak) . Again following the standard treatments of the problem we are going 
to neglect the barrier. For a discussion of its effect on the modes see (DeWitt 1975), section 5.2. Thus (40) becomes 

(42) 



9P 9r*2 

(in what follows the index in (^ is suppressed). Since V'o(^ = 0) must be finite it follows from (39) that 

<j){r = 0) = (43) 

Relations (42) and (43) define the problem. 
Relation (42) can be written in the form 

= (44) 



dudv 



Hence any function that depends only on m or u (or the sum of two such functions) is a solution of (44). We can 
now make contact with the accelerating mirror problem (DF,I). One set of modes satisfying (44) and the boundary 
condition (43) is given by 

(^^(u, v) = „ ^ — (exp(-iww) - exp {-iujp{u))) (45) 



Another set of modes satisfying the boundary condition is immediately obtained from (45) 

2 

i^a;(w, v) = — (exp {-iujf{v)) - exp {-iuju)) (46) 

The modes Lpuj{u,v) of (45) describe sinusoidal waves incident from r* = oo as it is clear from the sign of the 
exponential in the first term; the second term represents the outgoing part which has a rather complicated behaviour 
depending on the motion of the mirror. These modes constitute the in space and should obviously be absent before 
coUapse, 

a{uj') |0m) = (47) 

Similarly the modes ip^{u,v) describe sinusoidal outgoing waves travelling towards r* — od (emitted by the star) as 
can be seen from the exponential of the second term. Correspondingly the first term is complicated. 

Recall that apart from the modes (paj{u,v) we should include modes quj{u,v) that contain no outgoing component; 
instead they are confined inside the black hole (the letter q complies with (Hawking 1975, 1976)). The modes quj{u, v) 
are undetectable by an outside observer, but they are needed to make the set f>uj,quj complete. Correlations between 
the two sets is part of the information problem alluded to in the Introduction. 

The (poj {u, v) modes should be absent from the out vacuum 

a(uj) \Oout) = (48) 

The state \Oout) corresponds to the state where there are no outgoing particles detectable by an outside observer. 
The two representations are connected by the Bogolubov transformation 

/•oo 

a{uj)= duj{a{uj,uj')a{uj')+(3*(uj,uj')a'i{uj')) (49) 

Jo 

Using (49) and its hermitean conjugate we may immediately verify that the expectation value of the number of 
excitations of the mode ip^ in the \Oin) vacuum 



(Om| iV (w) |0m) = / dLj'\[3{Lj,Lj')f (50) 

Jo 

The matrix element (3{uj, lo') is given by (21) of I (we use z for either r or r* as the case may be to ease up on notation 
and make contact with I) 

(3{uj,uj') ^-i dzip^,{z,0)—ip^{z,0) + idz(—ip^,{z,0)jf>^{z,0) (51) 

The integration in (51) can be over any spaceHke hypersurface. Since collapse starts at t = the choice t — for the 
hypersurface is convenient. The in modes evaluated at i ~ are given by the simple expression (45) (i.e. p{u) — u) 

1 
ipuj (u, v) ~ — (cxp(— io^z) — exp (iuz)) 9 (z) (52) 



where the presence of (z) emphasizes the fact that z is a radial coordinate. The Cp modes are given by (46) with / 
depending on the history of the collapse. Relation (51) is rewritten in the form (the cndpoints of integration shall be 
stated presently) 

0(lo,uj'] = —i I di 

)'_ 
2V7iu'' 

The above expression may be rearranged in the form 

/3(w,w') = \= f " dz |e'"^'^ - e'''^'^ 6> (z) {Lo'e""^^' - ujfe-''^^] + ^'^ ~^ / dz ie'"^' ' - e-'"^' '] 9 {z)e''^' 

At:\/ujuj' Jo ^ -* ^ttvlulj' Jo ^ ^ 

(54) 

where the limits of integration are displayed (we took inro account the fact that the argument of / runs up to vh)- 
The rays with w < do not affect the amplitude due to the presence of the 9 function. The first and second integral 
in the above relation will be denoted by Pi{uj,lu') and /32(t^,w') respectively. Note that the second integration is of 
kinematic origin and totally independent of the collapse. We also quote the expression for the a{uj,Lu') amplitude 

a{uj,uj') ^ -i dzip^'{z,0)—ipl^{z,0) + i dz l—ip^,{z,0) j ip*^{z,0) (55) 

Observe the unitarity condition 

duj (a {uJi,uj) a* {102, uj) - (3 {lji,uj) (3* (w2, w)) = S {lji - W2) (56) 



P(lu,uj') = -i dz—^\e-'''''-e''^'']9(z)^^{f'e-''^f-e'''']+ (53) 

J 2V7rw' 1^ J 2^/1:10 

+i f dz "^^ L-"^'' - e'"^''-] 9 (z) —!^= \e-"^f - e'"^'] 

J 2a/W^ I J 2./^ ^ J 



which is a consequence of the fact that the set of in states is complete. 

Recall from I that we can introduce quantities A{uj,uj'), B{u},lu') that are analytic functions of the frequencies via 

, ,, A(lj,lj') ^. ,, B{uj,uj') , , 

aiu;,u;')^^^^,/3{ij,ij')^ "- "_' (57) 



The quantity B(uj,lv') is read off (53) (and A{u!,U!') from the corresponding expression for a(w,a;')). From the 
definitions of the Bogolubov coefficients, the explicit form (53) of the overlap integral and expressions (45) and (46) 
for the field modes one can deduce that 

B*{uj,uj') = A{-uj,uj'),A*{uj,Lu') == B{-uj,uj') (58) 

The above relations allow the calculation of a{uj,uj') once P{uj,uj') is determined. 



IV. CALCULATION OF THE BOGOLUBOV AMPLITUDES 

The strategy we adopt in handling (54) is as fohows. The first integral will be evaluated via an asymptotic expansion 
in negative powers of uj' ^ which will in fact show that the uj' integration in (3) is logarithmically divergent. We start 
with the second integral in (54): 






(59) 



The integrals in (59) are readily evaluated in terms of the function ( and its complex conjugate (* (see e.g. Heitler 
(1954), pages 66-71): 

f°° 1 

C(a;) ^ -i / e"'''dK = P-- iT:5{x) (60) 

Ja X 

Since we are interested in the asymptotic limit w' ^ oo the argument of the 5 function in (60) never vanishes so it is 
only the first term in (60) that is operative as far as the calculation of the /32(ci;, w') goes. (The 5 proportional term 
is relevant in the calculation of the a{Lo, lo') amplitude via relations (57), (58).) Thus asymptotically in the said limit 

P^{u;,J)^—^ (61) 

We turn to the first integral /3i(a>,ti>') in (54). Rather than splitting it to four integrals we perform an integration 
by parts to get (54) in the form 

1 /, ./ rVH 1 1 

/3i(u>,u>')- \— dze-*'"/(")-*" " + ^ sin (u> V) e-''"-^^"") (62) 

27r V u! Jo 27r Vww' 

This is the same integration by parts that was used to obtain (35) of I and also the one that is used in DF to go 
from their (2.10a) to (2.10b). The lower limit contribution of the integration by parts vanishes (being proportional to 
sinw'z), quite irrespectively of the value of /(O). The second term in (62) (originating from the upper limit) oscillates 
rapidly since the exponent tends to infinity. Thus the term tends distributionally to zero it may be neglected as in 
DF (also it is one power of lj' down compared to the first term) . So asymptotically in lo' we are entitled to write 

/3i(c.,c.') c. ^^\[^ r dz^^-f^^^-^-'^ (63) 

27r V ui Ja 

To bring the singularity in the integral to zero we make the change of variable 

z = vh — X (64) 

and rewrite j3i{ijj,uj') in the form 



/3i(u;,a;')-— ^^y- / dxe-'^^'^^^e''^'^ (65) 

where the function q{x) = f (vh — z) is defined in the range < x < vh and has the properties that follow from (37), 
(38): 

g{x):^-AM\n(-Yx^O (66) 

givH) - /(O) (67) 
We isolate the integral 

/ = / " (ixe-^'^»(^)e^'^'^ (68) 



To obtain the asymptotic behaviour of (68) for lo' large we adopt the standard technique of deforming the integration 
path to a contour in the complex plane; see (Bender and Orszag 1978), chapter 6; (Ablowitz and Fokas 1997), chapter 
6; (Morse and Feshbach 1953), p. 610 where a very similar contour is used in the study of the asymptotic expansion 
of the confluent hypergeometric. The deformed contour runs from up the imaginary axis till iT (we eventually take 
T ^ oo), then parallel to the real axis from iT to iT + w/f , and then down again parallel to the imaginary axis from 
iT + V}i to u/f . The contribution of the segment parallel to the real axis vanishes exponentially in the limit T ^ oo. 
We thus get 

/>CO r-OO 

I = i dse-'^'''e-"^<^('-'^ - i / dse"^' '^'^ "+'■''> e~"^<^('' "+'■''> (69) 



In both integrations the dominant contribution comes from the region where s ~ 0. The second integral (including 
the minus sign in front) takes the form using (67) 

/-oo ii^'vHp-iojf(0) 

^ie^'^^H^-i.^fia) dse-^' = -i (70) 

Jo ^ 

In the first integral we use the asymptotic form (66) valid for small x and write 

is\\ t is\ / n \ / s \ «4Mw 



expi-iujg{is)) = exp iiUdcjln ( _ 1 1 = ( _ 1 = exp (^--4Mwj (^-j 

where we took the branch cut of the function x^'^ to run from zero along the negative x axis, wrote x^' 
exp {iuj (In a; + i2N'K)) and chose the branch iV = 0. Thus 






f /IT \ f / TT \ 



^...,Az.-.4M. r(i + »4AM 

(71) 



We substitute (70) and (71) in (69) and then in (65) to get 






Collecting (61) and (72) we get 

^.-.4M^ r(l + z4Mu>) e— /W 1 



/^ (-,-') - -^^^^exp (-:^4M.) i^-^"" ^^::.:Z ^ - ^?^^=, + :r-^ (73) 



The day is saved by (37) which causes an exact cancellation of the last two terms. The first term on its own 
immediately leads to the black body spectrum. Taking its modulus, squaring, and using the property 

|r(l + iy)|^ = Try/ sinh (Try) 

we get 

,2 4M 1 



i/3(^,^')r 



27ra;' e 



I p%TXU>M 



V. CONCLUSION 

The objective of the paper was to prove that the Bogolubov amplitude P{lj,lj') has the asymptotic form (4) and 
that the radiation emitted has the spectrum of a black body. Standard quantum mechanics dictate that we should 
specify the initial and final states before calculating the transition amplitude, and to this end we considered the 
photon in states before the collapse and the out states after collapse has taken place (the final state does not have 
a qLj{u,v) component; cf the remark following (47)). We consider the case where the gravitational field is initially 
weak so that we may use the simple modes (52) as in modes; this is standard practice and leads to expression (65) 
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for the amplitude. Unfortunately the term (59) and the consequent (61) are missing from the standard treatments. 
Term (59) does not depend on the kinematics of the collapse and is certainly there in order to give a diagonal S 
matrix in the trivial case where collapse is never initiated and nothing is produced. A correct treatment of the large 
Lu' asymptotics of the amplitude (65) yields one contribution that leads to the black body spectrum and a second 
contribution (which again is missed in the standard treatments; see Appendix A) that precisely cancels (61). However 
this second contribution does depend on the collapse and the exact cancellation takes place only in the case of an 
initially weak gravitational field. In the light of the above remarks the black body result is indeed idependent of 
the details of the collapse (as often asserted) but does depend on the assumption of an initially weak gravitational 
field. The question as to what would happen in the case of an initial photon state in a gravitational background 
corresponding to an advanced state of collapse is not answered either here or in the standard treatments of Hawking 
radiation. The question may be of academic interest in the context of gravity, but it may be relevant in other cases 
where the analog of Hawking radiation is expected to occur. 

The derivation of the black body spectrum presented in this note is based on the calculation of the Bogolubov 
(3 {ll!,uj') amplitude. As emphasized in the Introduction this quantity is by definition time- independent, and thus the 
question as to where and when the photons are produced simply does not arise. It is certainly true that were it not 
for the singularity on the horizon the thermal spectrum would not arise. There are of course arguments based on 
calculation of local field quantities to the effect that the regime in question is the important one; or even simpler 
classical arguments related to the red shift during collapse ((Weinberg 1972), p. 347). However such statements may 
be misleading in connection to the quantum mechanical calculation of global (time independent) quantities. Similarly 
attempts to distinguish between "transient" and "steady state" radiation at the level of the a and /? amplitudes are 
bound to fail; the emphasis in the literature on the behaviour of the amplitude near the horizon has unfortunately 
led to such statements. One of the main conclusions of this note is that the correct derivation of the thermal result 
requires the consideration of the function /(f) throughout its range and not just of its asymptotic part. 
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Appendix A: Comparison with previous derivations 

The black body result is often obtained via a sequence of somewhat peculiar mathematical steps. One often starts 
(see e.g. (Birrell and Davies 1982) p. 108 and also references cited therein) with expression (65) thus wrongly 
disregarding (61). Following that step one unaccountably uses (66) throughout the range of integration and not just 
asymptotically near the horizon where it is valid. One thus (incorrectly) gets dropping constant prefactors 

(3i{^,cj')^J dx[-) e-- (74) 

That (74) is a wrong approximation to the original integral may easily be seen by the fact that the use of (66) has 
changed the behaviour of the integrand at a; = vh (the non-singular end) and an asymptotic estimate similar to that 
given in section 4 would not lead to the cancellation that took place between the last two terms in (73). One then 
proceeds to rewrite (74) by rescaling uj'x —* x 

Pi[uj,uj ) Ki h (LU ) / axx e (lb) 

Jo 

Since one is chasing the ultraviolet divergence one simply sets uj'vh — oo, changes variable p — ia and rotates in the 
complex plane to get (75) in the form 

/•oo 

Pi{uj,uj') w ^-^-4M (^,^-^^4M-l g-5^4M / ^ac-" a'^ (76) 

Note that setting lo'vh — oo certainly does not amount to a systematic expansion in {lo') . The a integration yields 
r(l + iiS) and one thus obtains the form for the /3i amplitude leading to the black body spectrum. On the other 
hand the step from (75) to (76) is again incorrect. Integral (74) can be performed exactly in terms of the confluent 
hypergeometric function and the asymptotic estimate for large w' may be examined afterwards. Indeed let us rescale 
the variable in (74) x -^ x/vh and rewrite 

Jo 

= „»-'4M+i^-,^4M^_^^ (^ _^ ^^4^^ 2 + iLuAM, iLu'vn) 
iuj + 1 
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where M is the confluent hypergeometric. We can now examine the asymptotic hmit of (77) for large to' . The 
asymptotic hmit of the confluent M{a^h,i\z\) for large values of \z\ is given by item 13.5.1 of (Abramowitz and 
Stegun 1972) (z = iuj'vh)- In the case h — a + 1 some simplifications occur and we get 



M{l + iuj,2 + iuj,i\z\) 



(l + zc.)e^l^l^+zr(2 + z^)^^ 



(78) 



(other terms are down by higher powers of 1/ \z\). The second term of the above relation combined with the prcf- 
actors in (77) does feature the F (1 + ioj) e^^ factor characteristic of the black body spectrum. The reason for the 
discrepancy between (76) and (78) lies in the fact that one should first evaluate the integral in terms of the confluent 



and then take the 



oo limit rather than take the limit first. The rotation in the complex plane stumbles upon 



the Stokes phenomenon for the confluent (different limits for \z\ -^ oo depending on argz)). In short the black body 
spectrum (76) is obtained by (a) incorrectly dropping (61), (b) incorrectly approximating (65) by (74), (c) wrongly 
estimating the oj' asymptotics of the latter. 

Appendix B: Technical remarks on the collapse of a sphere of dust 

We write down the expression for the Regge- Wheeler coordinate of the stellar radius as given by (11), (16): 

R* iv) = ^ (1 + cosr/) + 2Mln [^ (1 + cost?) - 1 
It takes a sequence of trivial trigonometric transformations to bring this to the form 



R* (t]) = - (1 + cosry) + 2Mln 



sin 22+2 sin mip. 
cos2f 



(79) 



We similarly transform the argument of the logarithm in (14): 

'1 



t{ri) = 2a\/l - ka? 



(rj + sinry) + ka rj 



2Mln 



no+ri 



sm ■ 



(80) 



Both R* {rj) and t (77) diverge as 77 ^ 770 (i.e. as the stellar radius approaches the Schwarzschild radius), but the 
combination v = i + R* does not: 



2ay/l - ka? 



— (77 + sin rj) -\- ka 77 



a cos^ ? + 2M In 



C0s2f 



The quantity C2 in (21) is the derivative of the above evaluated at tjq: 



C2 = 2a\/l - ka^ ( x + ;^ COS770 



■ sin 770 + 2M cot 770 



The quantity Ci in (21) is trivially obtained by differentiating (12) 

61 =: -sm77o 

To calculate A in (18) we rewrite the latter in the form 

R~2M 



2M 



j^^-i/2M 



We take logarithms of both sides of (84) to rewrite it in the form 

, (R-2M\ , ^ i 

We now use (9), (13), (85) and a sequence of trigonometric identities to obtain 



(81) 



(82) 



(83) 



(84) 



(85) 



In 



■ sm 7/0 



In tI Jl - ka^ 

M 



{rjo + sin 770) + ka^r/o 



(86) 
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thus determining A. The crucial step in the calculation of A lies in the fact that although both sides of (85) diverge 
in the limit 77 — > 770, there is precise cancellation of a term 

In sm 



on each side thus leading to the finite result (86). 

For a given value of 770 the corresponding value t{rio) given by (12). We take advantage of the fact that in the 
comoving frame the endpoints B and A of the ray BCA lie at r = a. Thus the unknown is the time tn in the 
comoving frame where ray 11 hits the star's surface. The latter may be determined via the equation of a null geodesic 
inside the star obtained through (5). The path of a light ray BCA propagating inside the star is given by 

*« dt r dr 

(87) 



To evaluate the left hand side of (87) we change variable according to (12) 

dr 1 . , 

(1 + cos 77) 



dr) 2Vfc 

Using (11) for R{t) the left hand side of (87) immediately yields 

1 



The right hand side of (87) yields 



('70 - Vh) (88) 



— = arcsin ( Vka 
Vk V 

which in the limit (10) reduces to 2a. Combining with (88) we get 

m-VH = sVka (89) 

In the limit (10) we get that vh lies very close to vq. 
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